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We show that driven dislocation assemblies exhibit a set of dynamical phases remarkably similar 
to those of driven systems with quenched disorder such as vortices in superconductors, magnetic 
domain walls, and charge density wave materials. These phases include jammed, fluctuating, and 
dynamically ordered states, and each produces distinct dislocation patterns as well as specific fea- 
tures in the noise fluctuations and transport properties. Our work suggests that many of the results 
established for systems with quenched disorder undergoing depinning transitions can be applied to 
dislocation systems, providing a new approach for understanding dislocation pattern formation and 
dynamics. 

PACS numbers: 74.20. Mn, 74.72.-h, 71.45.Lr, 74.50. +r 



There are numerous examples of systems of collectively 
interacting particles that, when driven externally, depin 
and undergo dynamical pattern formation and/or dy- 
namic phase transitions, such as a transition from a fluc- 
tuating to a nonfluctuating state. Such systems include 
domain wahs, driven vortices in tyjpe-II superconductors 
[l|-[6[ , sliding charge density waves \?\ , and driven Wigner 
crystals [8]. In these systems, fluctuating and intermit- 
tent dynamics arise just above depinning when an applied 
external force is increased from zero, while for higher 
drives the particles dynamically order into patterns such 
as anisotropic crystals or moving smectic phases with dif- 
ferent types of fluctuation statistics (9l-[l3|. Dislocations 
in materials are known to undergo a transition at the on- 
set of irreversibility or yielding that has similarities to de- 
pinning [14, 15]; however, it has not been shown whether 
driven dislocations can exhibit the same general features 
as other systems with depinning transitions. Establish- 
ing such a connection could potentially open an entirely 
new paradigm for understanding driven dislocations. 

It is known that organized dislocation structures 
within individual crystals, such as tangles, cells, or pla- 
nar walls, can become more refined and better defined 
as stress or strain increases. Two-dimensional (2D) and 
three-dimensional (3D) dislocation dynamics simulations 
based on linear elasticity theory predict self-organization 
of dislocation assemblies into varying configurations, such 
as pileups near the yielding or depinning transition \l4- 
[3 and 2D mobile walls or 3D slip bands Eliol 

under an external drive. Below a critical stress where dis- 
locations show no net motion, the system is considered 
jammed 0,0, while intermittent or strongly fluctuat- 
ing behavior with highly jerky or avalanche- like motion 
occurs above the critical stress [l^- Avalanche behav- 
ior with power-law velocity distributions is pro pos ed to 
be a signature of critical dynamics 17, lH, 23, 2^. No 
correlations between the transitions in patterning and 
the intensity of applied stress or strain have been es- 
tabhshed before now. Here we demonstrate that driven 



dislocation assemblies exhibit the same nonequilibrium 
phases observed for collectively interacting particle sys- 
tems exhibiting depinning, including pattern organiza- 
tion in the pinned state, a strongly fluctuating intermit- 
tent phase with a coexistence of pinned and moving par- 
ticles [H, H, [lol, [H, [26|, and at higher drive, when the ef- 
fectiveness of the substrate is reduced, a phase in which 
the dislocations organize into moving wall structures [H- 
0, 0, 0, [lll - [l3| . The onsets of these different dynamical 
regimes are correlated with pronounced changes in the 
transport curves hj, l5| , noise properties [3|, 22} , and spa- 
tial structures 0, S Isl • The onset of these phases can be 
observed via changes in the dislocation structure, mobil- 
ity, velocity distribution, and velocity noise. Our work 
implies that many of the established results obtained for 
driven vortex and other systems can be used to under- 
stand dislocation dynamics. 

We utilize a discrete dislocation dynamics model with 
periodic boundary conditions for a 2D cross section of 
a sample containing Nd = 480 straight edge disloca- 
tions that glide along parallel slip planes. This model 
was previously shown to capture the behavior observed 
in stressed anisotropic materials, particularly the inter- 
mittent flow near the onset of motion [13, [in, [in. An 
equal number of positive and negative moving disloca- 
tions are randomly placed in the sample and can move 
in the positive or negative x-direction depending on the 
sign of their Burgers vector b. At most one dislocation is 
allowed to reside on a plane, so in-plane pile-ups are pre- 
vented. The dislocations are also restricted from leaving 
their assigned glide plane. Rather than imposing an an- 
nihilation rule 15, |l3, IJs^, we enforce that two adjacent 
glide planes must be separated by at least Sy^ where Sy 
is on the order of the Burgers vector of the dislocations 
[2ll,[2i,[30|. 

The dislocations interact via a long-range anisotropic 
stress field that is attractive between two oppositely 
signed dislocations and repulsive for liked-signed pairs. 
We utilize a replicated image model to allow a large num- 
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ber of dislocations to be simulated efficiently over long 
Within the simulation volume, 



times [31|J. Within the simulation volume, all disloca- 
tions are subject to the stress fields of all surrounding 
dislocations regardless of their position. To best make 
the connection with particle systems, nucleation of fresh 
dislocations during loading is suppressed. Under an ex- 
ternal applied stress Text, dislocation i moves along x in 
its assigned plane according to an overdamped equation 
of motion given by r]^ = bi {YJj^i ^int(rj - r^) - Text) 
where Xi is the x coordinate of dislocation i at point 
= {xi^Hi) with Burgers vector value 6^, t] is the ef- 
fective friction, and T\nt(yj — ^i) is the long-range shear 
stress on dislocation i generated by dislocation j. The ex- 
ternal load on a dislocation is proportional to the stress, 
Fd = bText' For r = {x,y) = {xj ,yj)-{xi,yi), rint(rj-ri) 
for an edge dislocation with Burgers vector value b is 
rint(r) = bii[x{x^ - ^^)]/[27r(l - u){x^ + y^)^] where /i 
is the shear modulus and u is the Poisson's ratio. The 
length of the square simulation cell L is set to unity and 
the simulation volume remains fixed throughout loading. 
We normalize our units such that 6 = 1, = 1, and 
/i/27r(l — i^) = 1. The system is initially allowed to relax 
without an applied external drive. After relaxation, the 
external drive is applied quasi-statically, with sufficiently 
long waiting times between increments to avoid transient 
effects. We measure the average absolute value of the dis- 
location velocities {\v\) as a function of the stress. This 
is analogous to the voltage versus applied current curve 
for vortices in superconductors. 

In order to detect and characterize the dislocation con- 
tent and charge of the wall structures, we examine the 
distribution of dx = \xi — the x-axis separation be- 
tween two dislocations lying on adjacent planes. The 
fraction of dislocation pairs with dx < w is denoted as 
= nd^<:w/^tot^ where w is the pre- assigned maximum 
wall width, ritot is the total number of pairs in this sys- 
tem, and nd^<w is the number of pairs satisfying dx < w. 
Here we set w = 0.05, although other reasonable values, 
such as w = 0.02, give qualitatively similar results. To 
distinguish unipolar from dipolar walls, we discriminate 
between those pairs of like and unlike sign that lie within 



the critical wall width. The difference B = P- 
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-P^ 



between the fraction of pairs of like (unlike) sign 

{P-\ ) is directly related to the net Burgers vector around 

one dislocation within the wall width w. Thus, when a 
dipolar wall forms, B approaches zero since the number 
of unlike sign pairs is almost the same as the number of 

like sign pairs. When a unipolar wall forms, is zero 

and B is equal to 

As the system of randomly positioned dislocations is 
allowed to relax under zero applied stress, the disloca- 
tions reassemble into a locked configuration determined 
by the long-range stress fields they collectively produce. 
The relaxed arrangement of the dislocations shown in 
Figure 1(a) is disordered and contains no percolating 
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FIG. 1: Stress maps of the sample range from large negative 
(blue) to large positive (red) stress, (a) The initial dislocation 
positions at zero load, (b) Just before yielding, the disloca- 
tions are predominantly located at pile-ups to form a single 
bipolar wall, (c) Above yielding at Fd = 3.6, the wall breaks 
apart and the structure exhibits intermittent dynamics, (d) 
At Fd = 8.0 there is a dynamical ordering into polarized walls, 
each composed of dislocations with the same Burgers vector 
orientation. 



walls. The internal stresses generated by this spatially 
random arrangement are high and are distributed uni- 
formly across the volume. 

For loads < Fd < 2.0, the dislocation pattern slowly 
changes after each load increment but {\v\) goes to zero in 
the long time limit, indicating that the system is in the 



jammed phase below the critical yield [2l|, |22|. Under 
these low drives, any dislocation motion merely causes 
the dislocations to lock into another immobilized pattern. 
Figure 1(b) illustrates a typical locked dislocation config- 
uration for loads just below critical yield (i.e., Fd < Fc), 
where most dislocations have assembled into a dipolar 
wall comprised of a disordered arrangement of positive 
and negative dislocations that cannot move past one an- 
other. Compared to the initial state [Fig. 1(a)], the inter- 
nal stress [Fig. 1(b)] remains high but is more localized, 
and large stress concentrations appear in the vicinity of 
the wall. Such walls are analogous to the model of a "po- 
larized" wall 32, 33|, with dislocations of predominantly 



one sign on one side of the wall and the other sign on the 
other side. They are thought to be responsible for the ob- 
served hysteresis in unloading or the Bauschinger effect 
in subsequent reverse loadings 0, [111 . Observations of 
polarized walls have been reported in crystals deformed 
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to large strains j36l-l38l| . 

Just above yielding, the dipolar wall structure breaks 
down as shown in Fig. 1(c) and the system enters a 
state characterized by strong fluctuations in the dislo- 
cation positions. The dipolar walls repeatedly break up 
and reform, while the remaining wall fragments become 
smaller at higher drives and show continual change. The 
fluctuating state persists up to Fd = 5.0, when a new 
type of dynamic pattern appears where the dislocations 
form unipolar walls composed of only one type of dislo- 
cation, either negative or positive, as shown in Fig. 1(d). 
These walls can be identifled as disordered tilt walls. 
An ideal model of a tilt wall involves a periodic array 
of edge dislocations that accommodate a tilt misorienta- 
tion between two adjoining crystals. The development 
of low-misoriented tilt walls is suspected to be a pre- 
cursor for the eventual formation of subgrains in heavily 
deformed crystals [H, lio] . Most importantly, when the 
unipolar walls form, the internal stress decreases in ex- 
tent and intensity. The alternating positive and negative 
stress pattern that develops along the wall in Fig. 1(d) 
is consistent with the theoretical prediction for an infl- 
nite a rray of perfectly aligned, like-signed edge disloca- 
tions [4l| from linear elasticity theory. Thus, the high 
applied drive enables the dislocations to assemble into a 
low energy ordered structure, a result that is consistent 
with the theories of substructure development proposed 
by Kuhlman-Wilsdorf [33]. 

In Fig. 2 we show that the changes in the dislocation 
structure produce signatures in {\v\) versus for the 
system in Fig. 1. The upper curve in Fig. 2(b) shows 
the simple linear dependence of {\v\) on Fd expected for 
a single dislocation. For the interacting system, {\v\) is 
zero below yielding for 0.0 < Fd < 2.0, increases nonlin- 
early for 2.0 < Fd < 7.5, and then starts to become linear 
again for Fd > 7.5. Figure 2(a) characterizes the ordering 
dynamics as a function of Fd. Since the system is initial- 
ized in a random state containing no walls, ^ at 

Fd = 0^ but as the load increases, P+_ reaches a maxi- 
mum just below the yielding point as shown in Fig. 1(b) 

where a large dipolar wall forms. Above yielding, 

decreases in the fluctuating regime when the walls break 
up, and it gradually drops to zero in the high-driving 
region where the unipolar walls form. To identify the 
formation of the unipolar dislocation walls, we measure 

P-\-\-, , which rises for Fd > 5.0 in Fig. 2(a). Also shown 

in Fig. 2(a) is the net Burgers vector of the walls, indi- 
cating that for Fd > 5.0 the walls are indeed unipolar 
and contain either exclusively positive or negative dislo- 
cations. Samples with smaller numbers of dislocations 



show the same general features in P.^ , P_ 



and 



their difference P, as shown in Suppl. Fig. 1 [42]. 

The overall dynamics illustrated in Fig. 1 and Fig. 2 
are remarkably similar to those observed in driven sys- 
tems with quenched disorder. For example, for vortex 
matter as a function of external drive, there is a low 
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FIG. 2: (a) P^ — (blue squares), the fraction of dipolar walls, 
vs Fd has a peak just below yielding. P__^++ (black circles), 
the fraction of uni-polar walls, passes through a plateau when 
the polarized wall state forms. B (red triangles) is a measure 
of the net Burgers vector in the walls, (b) The average abso- 
lute value of the dislocation velocity (solid lower curve) 
vs Fd. The upper dashed curve shows for non-interacting 
dislocations. Visible in the lower curve is a yielding point, a 
nonlinear region corresponding to the disordered or fluctuat- 
ing regime, and a linear region at high drives when the system 
is dynamically ordered. Points a, b, c, and d indicate the Fd 
values illustrated in Fig. 1. 



drive pinned phase, a strongly fluctuating phase where 
the vortex lattice structure is disordered, and a highly 
driven phase where dynamical pattern formation occurs 
[3, 0, [ill . The corresponding vortex velocity- force curves 
also show the same features: the fluctuating phase is cor- 
related with a nonlinear region, while in the dynamically 
reordered phase the velocity depends linearly on the drive 

The dynamical phases in the vortex system have also 
been characterized by changes in the velocity noise fluc- 
tuations across different regimes. Just above depinning 
in the fluctuating regime, there is a strong 1//" noise 
signal [l|, m associated with a bimodal velocity distribu- 
tion that indicates a coexistence of pinned and moving 
vortices [H, [lo|, EH . The onset of dynamical ordering is 
accompanied by a drop in the noise power 5*0 as well 
as the appearance of narrow band noise features [H, H^l • 
For the dislocation system, in Fig. 3(a) we plot the in- 
stantaneous velocity distribution P(|v|) in the fluctuating 
phase at Fd = 3.2. We flnd a bimodal velocity distribu- 
tion that appears because a portion of the dislocations 
are immobilized in pileups while other dislocations have 
broken out of pileups and are mobile. A similar bimodal 
velocity distribution appears in the fluctuating phase for 
driven vortices and colloids. The corresponding power 
spectrum 5'(/) presented in Fig. 3(b) of the time series 
of the average dislocation velocity has a 1//^"^ signal in 
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FIG. 3: (a) at Fd = 3.2 is bimodal in the fluctuating 

phase, (b) The corresponding S{f) from the time series of 
the velocity has a 1//^"^ shape, (c) for Fd = 8.0 in the 

ordered phase shows a sharp peak, (d) The corresponding 
S{f) has a characteristic peak indicating narrow band noise. 



this regime, in good agreement with studies of driven vor- 
tices. At Fd = 8.0 [Fig. 3(c)] in the dynamically ordered 
phase, P(l'^l) has a single sharp peak and the correspond- 
ing S{f) in Fig. 3(d) has a narrow band feature with a 
characteristic frequency generated by the formation of 
an ordered structure of unipolar walls. In Fig. 4(a) the 
noise power for a fixed frequency reaches a peak in 
the middle of the fluctuating disordered phase and then 
decreases as the dynamically ordered phase of unipolar 
walls is approached. 

By conducting a series of simulations for varied dislo- 
cation densities, p = Njj/L'^^ and analyzing the order- 
ing dynamics, we construct the dynamic phase diagram 
shown in Fig. 4(b). The lower curve indicates the yield- 
ing transition from the low drive jammed or pinned phase 
of dipolar walls to the fluctuating disordered phase. The 
onset of the dynamically ordered phase is defined as the 
force at which the unipolar wall structures start to form, 
and is plotted in the upper curve. As p increases, the 
yielding point rises to higher F^ since the dislocations 
have a more difficult time breaking through the dipolar 
walls that form. The increase in yield threshold with in- 
creasing p remains robust when we perform simulations 
with different initial dislocation configurations. Fig. 4(b) 
shows that the onset of the high drive dynamically or- 
dered phase also increases in a similar fashion with in- 
creasing p. This phase diagram exhibits the same fea- 
tures observed for vortex systems as a function of pin- 
ning strength vs external drive, where both the criti- 
cal depinning force and the onset of the ordering rise to 
higher drives with increasing pinning strength [13]. For 
the dislocation system, increasing p is equivalent to in- 
creasing the pinning strength. In our simulation, the pin- 
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FIG. 4: (a) So vs Fd at a fixed frequency of/ = 5xl0~^ 
has a peak in the fluctuating phase and drops as the dynam- 
ically ordered phase is approached, (b) The dynamical phase 
diagram Fd vs 1/p, where p is the dislocation density. The 
lower curve (red circles) indicates the onset of yielding and the 
upper curve (black squares) is the onset of the dynamically 
induced ordered phase; the fluctuating phase falls between 
the two curves. Both the critical yielding and the dynamical 
ordering shift to higher drives as the p increases. Inset: The 
same curves plotted on a log-log scale. 



ning strength is directly related to the dislocation dipole 
break stress given by /i/ (87rdy(l-z/)) [iHI , where dy oc 1/p. 
Thus, as Fig. 4(b) shows, the pinning strength scales as p. 
This differs from the recent simulation results obtained 
in smaller systems (21, 22|. 

The type of ordered state that forms in the strong driv- 
ing regime varies depending on the rate at which the 
external load is applied. The ordered polarized walls il- 
lustrated in Fig. 1 form under continuous sweeps of the 
load. If the load is instead instantaneously set to a high 
value, we observe a transient disordered phase followed 
by the formation of multiple lower density unipolar walls, 
as illustrated in Suppl. Fig. 2 [42[, instead of the two 
unipolar walls shown in Fig. 1(d). As for the nature of 
the yielding transition, recent work on depinning systems 
has suggested that plastic depinning falls into the class 
of absorbing phase transitions, specifically directed per- 



colation [25|, |29|, suggesting that the onset of yielding 
for dislocation systems could also fall into the class of 
directed percolation. 

In summary, we have shown that driven dislocation 
assemblies exhibit the same nonequilibrium phases ob- 
served for systems of collectively interacting particles 
such as vortices in disordered superconductors. These 
include a jammed phase analogous to a pinned state, a 
fluctuating or disordered phase, and dynamically ordered 
or pattern forming states. All of the states are associated 
with transport signatures such as changes in the trans- 
port noise fluctuations as well as features in the disloca- 
tion velocity vs applied shear, in analogy with velocity- 
force curves. 
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sured dynamical reordering curves for samples with dif- 
ferent numbers of dislocations Nd- The results in Fig. 
2(a) are for a sample with Nd = 480. Figure SI shows 




FIG. 2: (Supplemental) The red (blue) crosses are the lo- 
cations of the dislocations with positive (negative) Burgers 
vectors in a sample with Nd — 480 that has been instanta- 
neously loaded with Fd = 10. The time progression of the 
formation of multiple unipolar walls is illustrated, a Initial 
dislocation configuration, b After 1 time unit, c After 3 time 
units, d After 5 time units, e After 20 time units, f After 
149 time units. 



FIG. 1: (Supplemental) P-\ (blue squares), the fraction of 

walls containing dislocations with both positive and negative 
Burgers vectors, P__,++ (black circles), the fraction of walls 
composed of a single polarity of burgers vectors, and B, the 
net Burgers vector in the walls, vs Fd for samples with dif- 
ferent numbers of dislocations Nd- a Nd =48. b Nd — 96. 
c Nd = 144. d Nd = 192. e Nd = 240. f Nd = 360. The 
same three regimes, pinned, fluctuating, and reordered, occur 
for all values of Nd, but the location of the regimes shifts to 
different values of Fd as Nd changes. 



corresponding plots of 



-,++' 



and P as a func- 



tion of Fd for samples with Nd = 48, 96, 144, 192, 240, 
and 360. The three regimes, pinned, fluctuating, and 
reordered, illustrated in Fig. 2(a) also appear in the 
smaller samples illustrated in Figure SI, but the tran- 
sitions between the regimes shift to lower values of Fd as 
Nd decreases. 



SUPPLEMENTARY INFORMATION FOR 
"DYNAMIC PHASES, PINNING, AND 
PATTERN FORMATION FOR DRIVEN 
DISLOCATION ASSEMBLIES" 

C. Zhou, C. Reichhardt, C.J. Olson Reichhardt, and 
I.J. Beyerlein 

Size effects on dynamical transitions: We mea- 



Instantaneous loading: In a sample that is instan- 
taneously subjected to a high load, we find a transient 
disordered phase followed by the formation of multiple 
lower density unipolar walls instead of the two unipolar 
walls shown in Fig. 1(g). For the sample with Nd = 480, 
Figure S2 illustrates the time progression of this process 
when we instantaneously increase to = 10. 



